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I. INTRODUCTION 



Although the notion of a density matrix is one of the fundamental concepts discussed in the elementary courses of 
quantum mechanics, the structure of the set A^at of all density matrices of size TV is not easy to characterize Q, |^ Q . 
The only exception is the case TV = 2, for which A^2 embedded in M.^ has an appealing form of the Bloch ball. Its 
boundary, 9A^2, consists of pure states and forms the Bloch sphere. For larger number of states the dimensionality 
of Mn grows quadratically with N, which makes its analysis involved. In particular, for N > 2 the set of pure states 
forms a 2N — 2 dimensional manifold, of measure zero in the iV^ — 2 dimensional boundary OMn- 

In this work we compute the volume of Mn with respect to the Hilbert-Schmidt (HS) measure. The HS measure 
is defined by the HS metric which is distinguished by the fact that it induces the flat, Euclidean geometry into the 
set of mixed states. The (hyper)area of the boundary of the space of the density matrices, OMn, is also computed, 
as well as the area of (hyper)edges of this set - the HS volume of the subspace of density matrices of an arbitrary 
rank k < N. In the special case of A: = 1 we obtain a well-known formula for the volume of the space of pure states, 
equivalent to the complex projective manifold CP^~^. 

A similar analysis is also performed for the set of real density matrices. To calculate the volume of the set of 
complex (real) mixed states we use the volume of the unitary (orthogonal) groups and the volume of the complex 
(real) flag manifolds - these results are described in the appendix. 

A motivation for such a study is twofold. On one hand, the complex structure of the set of mixed quantum states 
is interesting for itself. It is well-known that for TV > 2 the D = N"^ — 1 dimensional set A^at is neither a Z?-ball 
nor a polytope, but, how it looks like? More like a ball or more like a polytope? Instead of using techniques of 
differential geometry and computing the average curvature on the boundary of the set Mn, we compute the volume 
of its boundary and compare it with the volume of the _D — 1 sphere, which surrounds the ball of the same volume as 
A4n- Such a comparison shows us, to what extend the shape of the body of mixed quantum states differs from the 
ball, in a sense that more (hyper)area of the surface is needed to cover the same volume. 

A complementary information characterizing the structure of a given set is obtained by calculating the ratio between 
the area of its boundary and its volume. Among all _D-dimensional bodies of a flxed volume, such a ratio is smallest 
for the Z?-ball. Hence computing such a ratio for the ZJ-dimensional body of mixed quantum states we may compare 
it with similar ratios obtained for D-balls, ZJ-cubes and D-simplices. 

On the other hand our investigations might be useful in characterizing the absolute volume of the subset of mixed 
states distinguished by a certain attribute. For instance, if ^describes a composite system, one may ask, what is the 
volume of the set of separable (entangled) mixed states 0,13. Furthermore, assume we are given a concrete mixed 
quantum state g. It is natural to ask, wether g is in some sense typical, e.g. wether its von Neumann entropy is close 
to the average taken over the entire set A4j\i with respect to the HS measure. To compute such averages (see e.g. @) 
it is usefuU to know the volume of Al at and to make use of integrals developed for such a calculation. 
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II. GEOMETRY OF Mn WITH RESPECT TO THE HILBERT-SCHMIDT METRIC 

The set of mixed quantum states M n consists of Hermitian, positive matrices of size N, normaUzed by the trace 
condition 

Mn ■■= {g ■■ g= g''; g>0; tig=l; dim{p) = N}. (2.1) 

It is a compact convex set of dimensionahty D = N"^ — 1. Any density matrix may be diagonaUzed by a unitary 
rotation, 

g^UAU-\ (2.2) 

where A is a diagonal matrix of eigenvalues A^. Due to the trace condition they satisfy ~ 1' the space of 

spectra is isomorphic with a {N — l)-dimensional simplex Ajv_i. 

Let i? be a diagonal unitary matrix. Since g = U BAB'^U\ in the generic case of a non degenerate spectrum 
the unitary matrix U is determined up to N arbitrary phases entering B. To specify uniquely the unitary matrix 
of eigenvectors U it is thus sufficient to select a point on the coset space Fl^^^ :— U{N)/[U{1)]'^ , called complex 
Rag manifold. The generic density matrix is thus determined by (N — 1) parameters determining eigenvalues and 
N'^ — N parameters related to eigenvectors, which sum up to the dimensionality D of A^at. Although for degenerated 
spectra the dimension of the flag manifold decreases (see e.g. H,l3)i these cases of measure zero do not influence 
the estimation of the volume of the entire set of density matrices. Several different distances may be introduced into 
the set Mn, (see for instance 0,I3)- ^^^i^ work we shall use the Hilbert- Schmidt metric, which induces the flat 
geometry. 

The Hilbert-Schmidt distance between any two density operators is defined as the Hilbert-Schmidt (Frobenius) 
norm of their difference, 



DHs{gi,g2) = \\gi - £'2||hs = V'Tr[(6ii - £'2)^]. (2.3) 

The set of all mixed states of size two acquires under this metric the geometry of the Bloch ball B'^ embedded in 
R'^. Its boundary, c^B'^ — contains all pure states and is called Bloch sphere. To show this let us use the Bloch 
representation of a A^ = 2 density matrix 

g=^+T-X, (2.4) 

where A denotes the vector of three rescaled traceless Pauli matrices {cTx, o'y, <Jz}/V^- They are normalized according 
to trAf ~ 1. The three dimensional Bloch vector t is real due to Hermiticity of g. Positivity requires trg'^ < 1 
and this implies |t| < 1/a/2 =: R2. Demanding equality one distinguishes the set of all pure states, g-^ — g, which 
form the Bloch sphere of radius R2. Consider two arbitrary density matrices and express their difference gi — g2 in 
the representation H2.4|) . The entries of this difference consists of the differences between components of both Bloch 
vectors fi and T2. Therefore 

D]is{gTi,gT2) = DE{Ti,f2) , (2.5) 

where De 'iS the Euclidean distance between both Bloch vectors in M^. This proves that with respect to the HS metric 
the set A^2 possesses the geometry of a ball B'^ . The unitary rotations of a density matrix g ^ U gW correspond to 
the rotations of r in M^. This is due to the fact that the adjoint representation of SU{2) is isomorphic with 50(3). 

The Hilbert-Schmidt metric induces a flat geometry inside M n for arbitrary A^. Any state g may be represented by 
(|2.4|) , but now the A represents an operator valued vector which consists of D — N'^ — 1 traceless Hermitian generators 
of SU{N), which fulfill trA^Aj = Sij. This generalized Bloch representation of density matrices for arbitrary A^ was 
introduced by Hioe and Eberly 91, and recently used in ^^l- The case A = 3, related to the Gell-Mann matrices, is 
discussed in detail in the paper by Arvind et al. j^. The generalized Bloch vector r (also called coherence vector) 
is D dimensional. In the general case of an arbitrary TV the right hand side of H2.5|l denotes the Euclidean distance 
between two Bloch vectors in M.^ Positivity of p implies the bound for its length 



|f I < I?Hs {VN, I V-) (V'l) = y ^-Rn- (2.6) 

In contrast to the Bloch sphere, the complex projective space CP^~^ which contains all pure states, forms for 
N > 2 only a measure zero, simply connected 2(A — l)-dimensional subset of the A^^ — 2 dimensional sphere of radius 
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Rn embedded in Thus not every vector r of the maximal length i?jv represents a quantum state. This is 

related with the fact that for N > 3 the adjoint representation of SU{N) forms only a subset of SO{N^ — 1), (see e.g. 
[T^V Sufficient and necessary conditions for a Bloch vector to represent a pure state were given in for = 3, 
and in fl3] for an arbitrary TV. Furthermore, by far not all vectors of length shorter then Rn represent a quantum 
state, as not all the points inside a hyper-sphere belong to the simplex inscribed inside it. Necessary conditions for 
a Bloch vector to represent quantum mixed state were recently provided by Kimura 14]. On the other hand, there 
exists a smaller sphere inscribed inside the set A4n- Its radius reads 



R 



N 



where qn-i denotes any state with the spectrum {j^^, jt^jO)- 



(2.7) 



III. HILBERT-SCHMIDT MEASURE 



Any metric in the space of mixed quantum states generates a measure, inasmuch as one can assume that drawing 
random density matrices from each ball of a fixed radius is equally likely. The balls are understood with respect to a 
given metric. In this work we investigate the measure induced by the Hilbert-Schmidt distance ()2.3|) . The infinitesimal 
distance takes a particularly simple form 



(dsHs)' = Tr[(df?)2] 

valid for any dimension N. Making use of the diagonal form p = UMJ^^ we may write 
Thus (|3.1|l can be rewritten as 

N N 



(dsHs)' - + 2E(^'^ - A,f |(C/-M[/), 



(3.1) 
(3.2) 

(3.3) 



Since the density matrices are normalized, X]i3=i = 1; thus X^iLi "^^i = 0. Hence one may consider the variation 



N 



of the A^-th eigenvalue as a dependent one, dAjv = — X^iLi^ dA^, which implies 

N N-l N-l N-l 

Y,{dKf = + (E dAO' = E dA.g.,dA, 

i—1 i—1 i—1 — 1 



(3.4) 



The corresponding volume element gains a factor -y/detg, where g is the metric in the {N — 1) dimensional simplex 
Atv-i of eigenvalues. From H3.4fl one may read out the explicit form of the metric gij 



(3.5) 



It is easy to check that the spectrum of the A'^ — 1 dimensional matrix g consists of one eigenvalue equal to N and 
remaining N — 2 eigenvalues equal to unity, so that detg — N. Thus the Hilbert-Schmidt volume element is given by 



1— ' 


0" 




" 1 • 


• l' 





1 


+ 


1 • 


• 1 



dVns = VnY[ dA, [] (A, - Akf \ J] 2Re(C/-MC/),ani(f/~M[/)jfe| 

j=l j<k j<k 



and has the following product form 



dV = d^(Ai, A2, Atv) X dz/] 



Haar 



(3.6) 



(3.7) 



The first factor depends only on the eigenvalues A^, while the latter on the eigenvectors of g which compose the 
unitary matrix U. 
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Any unitary matrix may be considered as an element of the Hilbert-Schmidt space of operators with the scalar 
product {A\B) — TrA'^B. This suggests the following definition of an invariant metric of the unitary group U{N), 

N N N 

{dsf := -TT{U-'dUr - E \(U-'dUW = ^ |(C/-M{/),,f + 2 ^ \{U~'dUW . (3.8) 

jk=l j=l j<k=l 

This metric induces the unique Haar measure vuaar on U{N), invariant with respect to unitary transformations, 
^'Haar(W^) = i^HaaiiUW), whcrc W dcuotcs an arbitrary measurable subset of U{N). Integrating the volume element 
corresponding to H3.8|l over the unitary group we obtain the volume 

4.^(^,1), - (3.9) 
Integrating the volume element with the diagonal terms in (|3.8|) omitted (in that case the diagonal elements of U are 

,(^)^ _ Vol[[/(jV)] ^ (2^)^(^-i)/^ 
(27r)^ 1!2!---(A^~ !)!■ 



fixed by Uu > 0) we obtain the volume of the complex flag manifold, Fl^!^-' :— U{N)/[U{1)^] 



Both results are known in the literature for almost fifty years 15]. However, since many different conventions in 
defining the volume of the unitary group are in use 16, 17, 18, .19^ 20, 21, 22J we sketch a derivation of the above 
expressions in the appendix and provide a list of related results. 

Comparing formulae H3.6|l and 13. 8|) we recognize that the measure ly, responsible for the choice of eigenvectors 

of p, is the natural measure on the complex fiag manifold Fl^^^ = U{N)/[U{1)^] induced by the Haar measure on 
U{N). Since the trace is unitarily invariant, it follows directly from the definition (|3.1|l . that the volume element with 
respect to the HS measure is invariant with respect to the group of unitary rotations, dVHs(f?) = dVusiU gU'^). Such 
a property is characteristic tq_any product measure of the form H3.7|l . Several product measures with different choices 
of n were examined in ^ j2^ |23] . 

Integrating the volume element H3.6|) with respect to the eigenvectors of g distributed according to the Haar measure 
one obtains the probability distribution in the simplex of eigenvalues 

N N 

p(|)(Ai, . . . , Aat) = C«s<5(l - ^ A,) n(A, - Afe)^ (3.11) 

j=i j<:k 

where for future convenience we have decorated the symbol P with the superscript consistent with the exponent 
in the last factor. As discussed in the following section the normalization constant may be expressed in terms 
of the Euler Gamma function T{x) 

= ■ (3.12) 

The above joint probability distribution, derived by Hall defines the measure /ins in the space of diagonal 
matrices, and the Hilbert-Schmidt measure H3.6|) in the space of density matrices A4n- 

Interestingly, the very same measure may be generated by drawing random pure states |0) g Hi 1^0.2 of a composite 
N X N system according to the Fubini-Study measure on CP^ Then the density matrices of size N obtained by 
partial trace, g = tr2(|(/))(0|), are distributed according to the HS measure [EEEES- Alternatively, one may generate 
a random matrix A of the Ginibre ensemble, (non-Hermitian complex matrix with all entries independent Gaussian 
variables with zero mean and a fixed variance) and obtain a HS distributed random density matrix by a projection 
g = A"^ A/trA^ A |^. A similar approach was recently advocated by Tucci |2^J, who used the name 'uniform ensemble' 
just for ensemble of density matrices generated according to the HS measure. 

IV. VOLUME OF THE SET OF MIXED STATES 

For later convenience let us introduce generalized normalization constants 

N N 



POO 

M^--- dAi . . . dA^siY: A. - 1) n n - ^ 



f (4.1) 
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with a,P > 0. These constants may be calculated using the formula for the Laguerre ensemble, discussed in the book 
of Mehta El, 



f .A. . ^ c,A„ .M- g A.) n A?- n iA. - A/ - n r''"Z^;/i^-^"'' ] («> 

Substituting xf — Ai we may bring the latter integral to the Gaussian form. Expressing it in the spherical coordinates 
we get the integral (|4.1|) and eventually obtain 

1 1 -,^rr[l+j/3/2]r[«+(j-l)/3/2] 



C^^^m ■ r[a7V + /3iV(iV-l)/2] -LiL r[l + /3/2] 



J-J-L T\l + B/2] \' ^ ' 



By definition = C'w^^ special case of the above expression reduces to (|3.12|l . 

To obtain the Hilbert-Schmidt volume of the set of mixed states Ai n one has to integrate the volume element (|3.6() 
over eigenvalues and eigenvectors. By definition the first integral gives while the second is equal to the volume 

of the flag manifold. To make the diagonalization transformation (|2.2|) unique one has to restrict to a certain order 
of eigenvalues, say, Ai < A2 < • ■ ■ < A^r, (a generic density matrix is not degenerate), which corresponds to a choice 
of a certain Weyl chamber of the eigenvalue simplex A^^i. In other words, different permutations of the vector of 
N generically different eigenvalues A^ belong to the same unitary orbit. The number of different permutations (Weyl 
chambers) equals to A^!, so the volume reads 

y'i^ Vo1hs(>(^) = ^ y ^ (4.4) 

The square root stems from the volume element (|3.6f) . and the index refers to the general case of complex density 
matrices. Making use of H3.12|l and (|3.1l)|) we arrive at the final result ^| 

^ ^/iV(2.)^(^^^)/^ (4.5) 

Substituting A^ = 2 we are pleased to receive v!f^ = 7r\/2/3 - exactly the volume of the Bloch ball B'^ of radius 
R2 = l/\/2. This result may be also found in the notes by Caves 22], who also derived an explixit integral for the 
volume of the set of mixed states for arbitrary A^. 

The next result V^'^^ = 7r^/(840-\/3) allows us to characterize the difference between the set A^s C and the 
ball B*. The set of mixed states is inscribed into the sphere of radius i?3 = \/2/3 « 0.816, while the maximal ball 
contained inside has the radius r^, = R3/2 sa 0.408. Using Eq. I|6.1|l we find the radius of the 8-ball of volume V3 
is p3 « 0.519. The distance from the center of Ai^ to its boundary varies with the direction in from rs to R3, 
in contrast to the N = 2 case of Bloch ball, for which R2 = r2 ~ p2 = l/'s/Z. The average HS-distance from the 
center oi AA^ to its boundary is equal to pa. Similar calculations performed for A^ = 4 give the maximal radius 
i?4 = Y^3/4 « 0.866, the minimal radius = R^/i sa 0.289 and the 'mean' radius p4 w 0.428 which generates the 
ball B^^ of the same volume as P4. In general, let denote the radius of a ball B^^^^ of the same volume as the 
set Mn- 

The volume Vn tends to zero if A" ^ 00, but there is no reason to worry about it. The same is true for the volume 
of the A^-ball, see (|6.1|l . This is just a consequence of the choice of the units. We are comparing the volume of an 
object in with the volume of a hypercube of side one, and it is easy to understand, that the larger dimension, 
the smaller is the volume of the ball inscribed into it. 



V. AREA OF THE BOUNDARY OF THE SET OF MIXED STATES 



The boundary of the set of mixed states is far from being trivial. Formally it may be written as a solution of 
the equation detp = which contains all matrices of a lower rank. The boundarycJA^jv contains orbits of different 
dimensionality generated by spectra of different rank and degeneracy (see eg. |3,I3). Fortunately all of them are of 
measure zero besides the generic orbits created by unitary rotations of diagonal matrices with all eigenvalues different 
and one of them equal to zero; A = {0, A2 < A3 < . • . < Ajv}. Such spectra form the N — 2 dimensional simplex 
A7v_2, which contains {N — 1)! Weyl chambers - this is the number of possible permutations of elements of A which 
all belong to the same unitary orbit. 
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Hence the hyper-area of the boundary may be computed in a way analogous to (|4.4|l . 

:= Vo1hs(M.) ^ ^^3^ ^J) ^ - (5-1) 

The change of the parameter a in 1)4. l|l from 1 to 3 is due to the fact that by setting one component of an N 
dimensional vector to zero the corresponding Vandermonde determinant of size N leads to the determinant of size 
N — 1 for (3 = 1 and to the square of the determinant for (3 = 2. Applying (|3.12() and H3.10|l we obtain an explicit 
result 

- f27r^^(^-i)/2 r(i)---r(jv + i) 

For N = 2 we get S^"^ = 27r - just the area of the Bloch sphere of radius R2 = l/\/2. The area of the 7-dim 

boundary of M3 reads 5"^ ^ = V2n^/105. 

In an analogous way we may find the volume of edges, formed by the unitary orbits of the vector of eigenvalues 
with two zeros. More generally, states of rank N — n are unitarily similar to diagonal matrices with n eigenvalues 
vanishing, A = {0, . . . , 0, A„+i < A„+2 < • • • < Ajv}. These edges of order n are N'^ — n? — 1 dimensional, since the 
dimension of the set of such spectra is — n — 1, while the orbits have the structure of U{N)/\U{n) x (C/(l))^~"] 
and dimensionality A^'^ — — (A^ — n). Repeating the reasoning used to derive H5.1|) we obtain the volume of the 
hyperedges 

{N~n)l c^i+2„,2) Vol(^^4"V ^ ^ ^ 

(2) (2) 
Note that for n = this expression gives the volume of the set A^^v, for n = 1 the hyperarea S}^ of its boundary 

OMn, for n > 2 the area of the edges of rank N — n. In the extreme case of ri = A^ — 1 the above formula gives 

correctly the volume of the set of pure states, (the states of rank one), Vol(CP^^^) = (27r)^^^/r(Af), see appendix. 



VI. THE RATIO: AREA/VOLUME 

Certain information about the structure of a convex body may be extracted from the ratio 7 of the (hyper)area of 
its boundary to its volume. The smaller the coefficient 7 (with the diameter of the body kept fixed), the better the 
body investigated may be approximated by a ball, for which such a ratio is minimal. And conversely, the larger 7, 
the less the body resembles a ball, since more (hyper)area is needed to bound a given volume. 

To analyze simple examples let us recall the volume of the A^-dimensional unit ball C and the volume Sn 
of the unit A^-sphere C 



vol(S^-i) nf 1 /27re\ , , 

BN:=voliB^)= \^ ^ ^ ^^^ _^ ^__(_) , (6.1) 



N r(f + i) V2^\N 



where the Stirling expansion |2J| was used for large A^. For small A^ we obtain the well known expressions. Si = 27r, 
S2 = 47r, S3 = 27r^, 5*4 = 87r^/3 and B2 = tt, B3 = 47r/3, B4 = iv^ 12. If the spheres and balls have radius L then the 
scale factor has to be supplied. In odd dimensions the volume of the sphere simplifies, vo^S^''"^^) = 2t& j {k — 1)!. 
Since the boundary of a A^-ball is formed by a A^ — 1 sphere, = S^~^, the ratio 7 for a ball of radius L reads 



vol(aB^) _ A^ 
vol(B^) ' T 



Intuitively this ratio will be the smallest possible among all A^-dimensional sets of the same volume. Hence let us 
compare it with an analogous result for a hypercube Dat of side L and volume L^. The cube has 2^ corners and 2A^ 
faces, of area L^~^ each. We find the ratio 

which grows twice as fast as for A'-balls. Another comparison can be made with simplices A^v, generated by (A^ + 1) 
equally distant points in R^. The simplex A2 is a equilateral triangle, while A3 is a regular tetrahedron. The volume 
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of a simplex of side L reads vol(A7v) = [L^ \/ {N + l)/2^]/iV!. Since the boundary of Aat consists of 1 simplices 
Aat-i we obtain 



,^ , vol(aAAr) / 27V N(N + l) ,^ 
vol(AArj y N + 1 L 

In this case the ratio 7 grows quadraticahy with N , which reflects the fact that simphces do have much 'sharper' cor- 
ners, in contrast to the cubes, so more (hyper)area of the boundary is required to cover a given volume. Furthermore, 
if one defines a hyper-diamond as two simplices glued along one face, its volume is twice the volume of A^r while its 
boundary consists of 2A'' simplices Ajv_i, so the coefficient 7 grows exactly as N"^ . 

Interestingly, the ratio 7 of the A^-cube is the same as for the A'^-ball inscribed in, which has much smaller volume. 
The same property is characteristic for the A'^-simplex. Hence another possibility to characterize the shape of any 
convex body F is to compute the ratio xi '■— vol[Bi(F)]/vol(F), and X2 '■= vol(F)/vol[B2(-F)], where Bi(F) is the 
largest ball inscribed in F while B2(F) is the smallest ball in which F may be inscribed. As stated above for cubes 
and simplices one has j{F) = 7[Bi(F)]. 

Such quotients may be computed for the rather complicated convex body of mixed quantum states analyzed with 
respect to the Hilbert-Schmidt measure. Using expressions (|4.4|l and (|5.1() we find 



7JV := \ . . V - = ^N{N-1) (N^ - 1). (6.5) 

The first coefficients read 72 = 3\/2, 73 = 8\/6, and 74 = 15vT2, so they grow with N faster than N^. A direct 
comparison with the results received for balls or cubes would be unfair, since here N does not denote the dimension 
of the set A^jv C M^. Substituting the right dimension, D — N"^ — 1, we see that the area/volume ratio for the 
mixed states increases with the dimensionality as 7 D^^^. The linear scaling factor L, equal to the radius Rn tends 
asymptotically to unity and does not influence this behaviour. 

Note that the set of mixed states is convex and is inscribed into the sphere of radius Rn, so for each flnite N the 
ratio 77V remains finite. On the other hand, the fact that this coefficient increases with the dimension D much faster 
than for balls or cubes, sheds some light into the intricate structure of the set tWjv. It touches the hypersphere 
of radius Rn along the 2N — 2 dimensional manifold of pure states. However, to be characterized by such a value of 
the coefficient 7 it is a rather 'thin' set, and a lot of hyper-area of the boundary is used to encompass its volume. In 
fact, for any mixed state £i G A^at its distance to the boundary OMn does not exceed the radius rjy ~ 1/iV. Another 
comparison can be made with the Z?-ball of radius L — tn = [N{N — 1)]"-*^/^, inscribed into A4n- Although its 
volume is much smaller than this of the larger set of mixed states, its area to volume ratio, 7 = _D/L is exactly equal 
to (|6.5(l characterizing /An- In other words, for any dimensionality N the set of mixed quantum states belongs to 
the class of bodies for which 7(F) = 7(61 (F)) holds. 

Using the notion of the effective radius pN, introduced in section IV, we may express the coefficients Xi for the 
set Mn SiS a. ratio between radii raised to the power equal to the dimensionality, D — N"^ — 1. The exact values of 
Xi = {i"n/pn)^ and X2 = {pn / Rn)^ , as well as their product x = X1X2, may be readily obtained from (|4.5|l . Let 
us only note the large N behaviour, x{-^n) — — 1)^^ +^ so it grows with the dimensionality D as D^^/^ while 
X(B^) = 1, x{^n) = N~^'\ and x(A,v) « N'^ . 

VII. REBITS: REAL DENSITY MATRICES 

Even though from physical point of view one should in general consider the entire set jv of complex density 
matrices, we propose now to discuss its proper subset: the set of real density matrices. This set, denoted by 
is of smaller dimension Di — N{N -|-l)/2 — 1 < D = N"^ — 1, and any reduction of dimensionality simplifies the 
investigations. While complex density matrices of size two are known as qubits, the real density matrices are sometimes 
called rebits ^2§\. In the sense of the HS metric the space of rebits forms the full circle B^, which may be obtained as 
a slice of the Bloch ball B^ along a plane containing 1/2. 

To find the volume of the set A4^ we will repeat the steps ()3.2|l - ()4.5|l for real symmetric density matrices which 
may be diagonalized by an orthogonal rotation, g — OAO^ . The expressions 

d£- = 0[dA + 0-MOA- AO-MOlO-i (7.1) 

and 

N N 

(dsHs)' =^(dA,)' + 2^(A,-A,)2|(0-idO),,f (7.2) 

1=1 i<j 
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allow us to obtain the HS volume element, analogous to (|3.6|l . 

Af~l 1...N 1---N 

d^Hs - n n iAj- - Afei • I n ^(o-mo),,i (7.3) 
j=i j<fc j<k 

As in the complex case the measure has the product form, and the last factor is the volume element of the orthogonal 

group (see appendix). Orthogonal orbits of a nondegenerate diagonal matrix form real flag manifolds Fl^^ — 
0(7V)/[C'(1)]^ of the volume 

voHF^M J - — ^ — - r[i/2]...r[7v/2] " ^^-^^ 

Here 0(1) is the reflection group Z2 with volume 2. 

The volume element ()7.3|l leads to the following probability measure in the simplex of eigenvalues 

N N 

P«(Ai, . . . , A^) = d^-'^S{l - 5] A,) n I A, -Ak\, (7.5) 

j=l j<k 

with the normalization constant given in (|4.3|) . Note the linear dependence on the differences of eigenvalues, in 
contrast to the quadratic form present in l|3.11|l . Taking into account the number A^! of different permutations of the 
elements of the spectrum A we obtain the expression for the volume of the set of Al* , 



which gives 



Vi,^> := VolHs {Ml) = ^ '7;:; ^ , (7.6) 



^(1) _ Vn 2^(2^)^(^-i)/4 r[^ ] ^^^^ 
m r[M-±i)] r[i] 



(1) _ viv . y.nj __Li^Zi nrfl + -l f77) 



As in the complex case we find the volume of the boundary of and in general, the volume of edges of order n 
with 0<n<A — 1. In the case of real density matrices these edges are Af(A^+ 1)/2 — 1 — n(n + 1)/2 dimensional, since 
the dimension of the set of such spectra is — n — 1, and the orbits have the structure of 0{N) /[0{n) x (0(1))^~"] 
and dimensionality A^(A^ — l)/2 — n{n — l)/2. In analogy to H5.3|l we obtain 



[N-ny. Vo1(p4"V " 

which for n — 1 gives the volume 5* of the boundary dAi^, and allows us to compute the ratio area to volume, 

voi(aA^«) _ Ni^T^r^ c^'^') 



The product of the last two factors is equal to the dimensionality of the set of real density matrices, Di = 
N{N + l)/2 — 1. Therefore, just as in the complex case, the ratio area to volume for Ai^ coincides with such a 
ratio 7 = Di/L for the maximal ball of radius L = tn = [N{N — 1)]~^/^ contained in this set. In the simplest case 
of Af = 2 we receive Vj'^'' = 7r/2 - the volume of the circle of radius R2 = l/\/2. The volume of the boundary, 
S — 7r\/2, equals to the circumference of the circle of radius R2 = l/\/2, and gives 7 = 2\/2 in agreement with (|7.9|l . 

VIII. CONCLUDING REMARKS 

We have found the volume V and the surface area S of the D — N"^ — 1 dimensional set of mixed states A4 n acting 
in the A-dimensional Hilbert space, and its subset containing real symmetric matrices. Although the volume of 
the unitary (orthogonal) group depends on the definition used, as discussed in the appendix, the volume of the set of 



9 



mixed state has a well specified, unambiguous meaning. For instance, for N = 2 the volume V2 may be interpreted 
as the ratio of the volume of the Bloch ball (of radius R2 fixed by the Hilbert-Schmidt metric) , to the cube spanned 
by three orthonormal vectors of the HS space: the rescaled Pauli matrices, {ax,o'y,o'z}/V^- 

On one hand, these explicit results may be applied for estimation of the volume of the set of entangled states 
H 0, m m 113, or yet another subset of M N. It is also likely to expect that some integrals obtained in this work 
will be usefull is such investigations. 

On the other hand, outcomes of this paper advance our understanding of the properties of the set of mixed quantum 
states. The ratio of the hyperarea of the boundary of Z?-balls to their volume grows linearly with the dimension D. 
The same ratio for ZJ-simplices behaves as , while for the sets of complex and real density matrices it grows with 
the dimensionality D as Z3'^/^ . Hence these geometrical properties of the convex body of mixed states are somewhere 
in between the properties of ZJ-balls and ZJ-simplices. 

Furthermore, we have shown that for any the sets of complex (real) density matrices belong to the family of sets, 
for which the ratio area to volume is equal to such a ratio computed for the maximal ball inscribed into this set. 

It is necessary to emphasize, that a similar problem of estimating the volume of the set of mixed states could be 
also considered with respect to other probability measures. In particular, analogous results presented by us in [3^ for 
the measure |l26i i3j| related to the Bures distance [ssl l36j allow us to investigate similarities and differences between 
the geometry of mixed states induced by different metrics. 

It is a pleasure to thank M. Kus for helpful discussions and to P. Slater for valuable correspondence. Finan- 
cial support by Komitet Badah Naukowych in Warsaw under the grant 2P03B-072 19 and the Sonderforschungsbe- 
reich/Transregio 12 der Deutschcn Forschungsgemeinschaft is gratefully acknowledged. 

APPENDIX A: VOLUMES OF THE UNITARY GROUPS AND FLAG MANIFOLDS 



Although the volume of the unitary (orthogonal) group and the complex (real) flag manifold, we use in our calcu- 
lations, were computed by Hua many years ago |l5| . one may find in more recent literature related results, which in 
some cases seem to be contradicting. However, different authors used different definitions of the volume of unitary 
group 0,0,0,0,113,122, so we review in this appendix three most common definitions and compare the results. 



1. Unitary group U{N) 

We shall recall (|3.8|l the metric of the unitary group U{N) induced by the Hilbert-Schmidt scalar product and used 
by Hua [H 

N N 

(d.)2:=-Tr([/-MC/)2=^|([/-MC/),,f + 2 ^ |(C/-M[/),,p , (Al) 

i=i j<fe=i 

which is left- and right-invariant under unitary transformations. The volume element is then given by the product of 
independent differentials times the square root of the determinant of the metric tensor. One has still the freedom of 
an overall scale factor for IjAljl which appears then correspondingly in the volume element. To keep invariance the 
ratio of the prefactors Cdiag and CoS of the diagonal and off-diagonal terms has to be fixed. Nevertheless one may 
introduce different scalings of the volume elements which we call dvA,dvB,dvc'- 

N 1---N 

di^A := |n(^"'dC/).. n ^ Re(C/-Mf/),feIm([/-Mf/),fc|; Cdiag = 1, = 2; (A2) 

i—1 j<k 

di/B := 2-^(^-1)/^ di^A; Cdiag = 1, Coff = l; (A3) 



dz/c := 2-^/2 dj.^. Cdiag = 1/2 Coff = l. (A4) 

The product in \A2\ . consistent with (|A1|) . has to be understood in the sense of alternating external multiplication 
of differential forms. Only the first convention l|A2p labeled by the index a was used in the main part of this work. 
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Note that the normahsation (|A4|) corresponds to the rescaled hne element (ds)^ = —^Tr(U ^dU)^. In general we 
may scale 



d^. := (cdiag)^/^(coff)^(^-^)/^dz.£ 



(A5) 



where the label x denotes a certain choice of the prefactors Cdiag and Cofi for diagonal or off-diagonal elements in 
(jAip . All these volumes correspond to the Haar measure which is unique up to an overall constant scale factor. Thus 
we deduce that 



VoU [UiN)] = 2^(^-i)/2Vo1b [UiN)] and Vole [UiN)] = [UiN)] . 



(A6) 



In order to determine the volume of the unitary group let us recall the fiber bundle structure U{N — 1) U{N) 



j2Af-l 



see e.g. |33. This topological fact implies a relation between the volume of the unit sphere S 



2N-1 



and the 



volume of the unitary group defined by the measure di's l|A3(l . for which all components of the vector {U ^dU)jk 
have unit prefactors, 



Yo\b[U{N)] = YoIb[U{N - 1)] X Vol[S2^-i]. 
To prove this equality by a direct calculation it is convenient to parametrize a unitary matrix of size N as 

-/it 



Un = 



J4> 







Un- 



h y/l-h<g) 



(A7) 



(A8) 



where (j) € [0,27r) is an arbitrary phase and /i is a complex vector with — 1 components such that \h\ < 1. This 
representation shows (we may arrange the two matrices in (jA8|l also in the opposite order) the relation 



U{N)/[U{1) X U{N - 1)] = CP 



(A9) 



since the second factor represents the complex projective space CP^ ^. In fact, if one calculates the metric (|A1I) we 
find 



(AlO) 



where (dsA?)^ means the metric for U{N) (the sign = shall indicate that we have omitted some shifts in (dsi)^ and 
(dsAT-i)^ that are not relevant for the volume) and (ds/J^ means the metric of the complex projective space CP^^^ 
with radius 1: 



(dsh)^ = dh^dh - 



(/itd/i + d/it/i)2 {h'fdh-dh'fh)^ 



4(l-|/ip 



(All) 



It is easy to see by diagonalizing this metric (eigenvalues 1 — 1/(1— |/ip), and otherwise 1) that the corresponding 
volume is that of the real ball B^^^^ with radius 1 and dimension 2N — 2. Thus one obtains 

Yo\x[U{N)] = Yo\x[U{N - 1)] x Volx[;7(l)] x c^-iVolp^^-^j ^ ^^^2) 

which for the measure ljA3|l with Cdiag = Coff = 1 reduces to (|A7|) . Applying this relation — 1 times we obtain 

Yo\b[U{N)] = Vol[S2^-i] X • • • X Vol[S3] X Vol[Si]. (A13) 

Taking into account that Vol[S^^~^] = /{N — 1)! and making use of the relation (|A6() we may write an explicit 
result for the volumes calculated with respect to different definitions (|A2l - IX4|l 



Volx[;7(iV)] = a 



X 



0!l!---(iV - 1)! 



(A14) 



where the proportionality constants read = 2^(^ i)/^^ 



-.u — 



1 and ag = 2"^/^. The result for Yo\a[U(N)] was 



rigorously derived in 15], Yo\b\U{N)] was given in |18|, while Yq\a\U{N)] and Yo\c\U(N)] were compared in [25 
In particular, VoU[C/(l)] = Vo1b[[/(1)] = 27r, while Volc[?7(l)] = \/2t: and Yq\a[U{2)] = Stt^, Yo\b[U{2)] = Att^, 
Yo\c[U{2)] = 27r3. 
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In general, the volume of a coset space may be expressed as a ratio of the volumes. Consider for instance the manifold 
of all pure states of dimensionality N. It forms the complex projective space CP^~^ — U{N)/[U{N — 1) x ?7(1)]. 
Therefore 



VolxiLP J - VoU[C/(l)] Vol^[C/(7V - 1)] ' ^^'^^ 



which gives the general result 



Yolx [CP^] =af^^ VoliB^'^] (A16) 



The scale factors read a^^ — 2^ and — = 1. For instance Vol^[CP"'^] = 27r which corresponds to the 
circle of radius \/2 , while Vo1b[CP^] = Volc[CP^] = tt, equal to the area of the circle of radius 1. The latter 
convention is natural if one uses the Fubini-Study metric in the space of pure states, Z?fs(|<p), IV")) = arccos(\/K), 
where the transition probability is given by k = Then the largest possible distance Dps — obtained 

for any orthogonal states, sets the geodesic length of the complex projective space to tt which corresponds to the 
geodesic distance of two opposite points on the unit circle, being identified. It is worth to add that Volc[CP'^] = 
Vol[S^*''+^]/Vol[S^] and this relation was used in 20] to define the volume Vole of complex projective spaces. We see 
therefore that different conventions adopted in ljA2l - IA4|I lead to different sizes (geodesic lengths) of the manifolds 
analyzed. 

Unitary orbits of a generic mixed state with a non-degenerate spectrum have structure of a (iV^ — A'^)-dimensional 
complex flag manifold Fl^^^'^ = U{N)/[U{1)]^ . Hence its volume reads 



{Wo\x[U{l)])'' l!2!..-(iV-l)i 



Volx[^r^] = / .^"^^ V = 4' ... (A17) 



with convention dependent scale constants a^' = 2^(^ i)/^ and aFg — = 1 [2^. It is easy to check that the 
relation 

Vo1x[f4^^] = Volx[CPi] X Volx[Cp2] X • • • X Vol^pP^^^] (A18) 



holds for any definition (|A2l - lA4)l . since the scale constants do cancel. 

For completeness we discuss also the group SU{N), the volume of which is not equal to Vol[C/(A^)]/Vol[C/(l)] 
0,0,123. To show this let us parametrize a matrix Y/v & SU{N) v 



Y, 



N 



e'-P 

e-^[^/(^-i)irjv-i 



h y/I-h(^h^ 



VW , (A19) 



where (/> e [0, 27r) is an arbitrary phase and his a, complex vector with — 1 components such that \h\ < 1. Condition 
detYAT = 1 implies TrY^^dY/v — 0. For instance, the metric IjAip gives, if the volume is concerned 

{dsf = -Tr(y-Ml/)2 - TiiW'^dWf. (A20) 

Since the first factor V is block diagonal the first term is equal to (d0)^A^/(A^ — 1) — Tr(Y^j';j^dY/v_i)^, while the 
second one gives the metric on CP^^^. Integrating an analogous expression in the general case of an arbitrary metric 
and using IjAlSp we obtain the following result 



^"HSUm - J:'rJ(T-''l)| iT&i Vol,[Sl/(iV- 1)1, (A21) 

which iterated — 1 times gives the correct relation 

with the stretching factor \/N. For instance, working with the measure (|A4p and making use of ljA14(l we obtain 
Yo\c[SU{N)] = y]v2(^-i)/27r(^+2)(^-i)/V[l!...(Af- 1)!], so in particular. Vole [5"!/ (2)] = 2tt^, Yo\c[SU{3)] = 
VStt^ and Vole [S'?7(4)] = V^tt^/S in consistence with results obtained in [T ll20ll2l l3i|. 
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2. Orthogonal group 0{N) 

The analysis of the orthogonal group is simpler, since (O^^dO)^ = —{O^^dO)^, so the diagonal elements of 
Tr(0~^dO)^ vanish. Thus we shall consider only two metrics (analogous to the measures (|A2I - IX4|) ') with different 
scalings, 

N 

(ds^)' -Tr(0-M0f = 2 ^ \{0-^dO)jkf, (A23) 

j<k=l 

used in section VII of this work, and 

(d.sf = (d.cf -^Tr(0-MOf ^ E 1(0-^0^.^, (A24) 

j<k=i 

which both lead to the Haar measure on the orthogonal group. 

To obtain the volume of 0{N) we proceed as in the unitary case and parametrize an orthogonal matrix of size N 

as 



O 



N 



Oi 

Oiv-1 



2 



h Vl-h(g)hT 



(A25) 



where Oi G 0(1) = ±, while h is here a real vector with iV — 1 components such that \h\ < 1. Representing 
the metric (dss)^ by these two matrices we see that the term containing only the vector h gives the metric of 
a real projective space. Integrating the resulting volume element (with scale factor 1) we obtain the volume of 
RP^~^, equal to iVol[S^~^]. Taking into account a factor of two resulting from 0(1) we arrive at VoIb[0{N)] = 
Yo\b[0{N — l)]Vol[S^^^], which applied recursively leads to 

VoIb[0{N)] = Vol[S^-i] X ... X Vol[Si] x Vol[SO] ^ H (^26) 

where Vol[S°] = Vol[0(l)] = 2. To get an equivalent result for the metric ljA23|l we have to take into account the 
factor \/2 which occures for each of N{N — l)/2 off-diagonal elements. Doing so we obtain 



N 1,^2 

VoU [0{N)] = 2^(^-i)/4Vob [0{N)] = 2^(^+3)/4 TT , (A27) 



in consistency with Hua 15]. In particular, VoU[0(l)] = Vo1b[0(1)] = 2, while VoU[0(2)] = AV2tt, Vo1s[0(2)] = 4tt, 
and VoU[0(3)] = 32V27r2, Vols[0(3)] = IGtt^. 

In full analogy to the unitary case we obtain the volume of the real projective manifold 

"^^'^ ^ - Vol.[0(l)]Vol,[0(iV-l)] • 

For the metric ljA24p this expression reduces to VolspP*^] = iVol[S'^]. Hence, this metric may be called 'unit sphere' 
metric, while the convention (jA23|) may be called 'unit trace' metric. 

In the similar way we find the volume of the real flag manifolds, used in analysis of real density matrices, 

(Volx[0(l)]) ^ 
Exactly as in the complex case we observe that the relation 

Volx [f4^^] = Volx [KP^] X Vo\x [MP2] x . . . x Vol^ [MP^-^] (A30) 



is satisfied for any definition of the metric. 
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Computation of the volume of the special orthogonal group SO{N) is much easier than in the complex case, since 
there are no diagonal elements in the metric and hence, no stretching factors. For any normalization one gets 

YolxlSOiN)] = ^yl^^^f^H^I = \ Volx[0(A^)] (A31) 

In particular, we get Vo1b[5'0(2)] = 2tt and Vo1b[5'0(3)] = Vole [5*0(3)] = Stt^. The latter results seem to be 
inconsistent with Volc[5'C/(2)] = 27r^, since there exists a one to two relation between both groups. This paradox 
is resolved by analyzing the scale effects jl^: the volume of SU{2) is two times larger than the volume of the real 
projective manifold conjugated to SO{i) of the appropriate geodesic length, VolcpP"^] = tt^. 
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